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convergence of the CFIE with multiscale meshes 
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Abstract — A hierarchical quasi-Helmholtz decomposition, orig¬ 
inally developed to address the low-frequency and dense- 
discretization breakdowns for the EFIE, is applied together with 
an algebraic preconditioner to improve the convergence of the 
CFIE in multiscale problems. The effectiveness of the proposed 
method is studied first on some simple examples; next, test on 
real-life cases up to several hundreds wavelengths show its good 
performance. 

Index Terms —Multilevel systems. Moment methods. Hierar¬ 
chical systems, CFIE 


1. Introduction 

The Electric Field Integral Equation (EFIE) is successfully 
used to numerically simulate the electromagnetic behaviour 
of structures in the frequency domain. Its Galerkin discretiza¬ 
tion using the classical Rao-Wilton-Glisson (RWG) functions 
defined in Q, however, produces an ill-conditioned matrix 
when the frequency is low or the discretization density is 
high. A hierarchical quasi-Helmholtz decomposition can help 
overcoming this issue, recasting the linear system in a form 
that can be easily preconditioned with a diagonal rescaling 
and an algebraic preconditioner on the fewer functions at 
the coarsest level. The EFIE suffers from other convergence 
problems, mainly related to resonances in closed bodies. A 
possible cure is to use the Combined Field Integral Equation 
(CFIE) instead, which mixes the EFIE with the Magnetic Field 
Integral Equation (MFIE), and exploits the fact that resonances 
occur at different frequencies for the two formulations (e.g. 
0). One is then tempted to use both the approaches and 
to apply the hierarchical quasi-Helmholtz decomposition to 
the CFIE, hoping to gain both the improvements. Focusing 
on the mesh density issue, we show that this can be done 
indeed, with simple examples and realistic cases to illustrate 
its effectiveness. Preliminary investigations along these lines 
were presented in the conference paper 0- 

Finally, we observe that the work in 0 shares our 
endeavour, yet with a completely different approach. 


H. Background 

The scattered electric and magnetic fields due to a time- 
harmonic surface current distribution on dVt, the boundary of 
the region in a medium with permittivity e and permeability 
/i are given by 
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respectively, where the neglected time factor is exp(jco’t) 
and g{r,r') = exp(—j/c|r — r'|)/47r|r — r'\ is the scalar 
Green’s function relevant to the surrounding medium, whose 
wavenumber is k. When considering a perfect electric conduc¬ 
tor with impinging time-harmonic electric and magnetic fields 
e* and h^, the unknown induced current j can be determined 
writing an equation on the boundary dVt, the EFIE 


— n(r) X n(r) x (e^(r) + e*(r)) = 0,r G dO.. (3) 


If the region is closed, the unknown induced current j can de 
determined writing an equation as approaching the boundary 
90 from the exterior, the MFIE 

— n{r) X {h^{r) + h\r)) = j{r), r -a- 90. (4) 

Expanding the unknown current in a finite basis, as the RWG 
defined in Q, and testing 0 and 0 against the same 
functions, in a Galerkin or Method of Moments (MoM) fashion 
(e.g. 0), we obtain two systems of linear equations 

[Z][j] = [e% [B][j] = [h^] (5) 
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Solving either one of 0 with a direct or iterative method 
gives the searched coefficients. To have better convergence 
properties, especially when considering closed objects, a linear 
combination of the two systems in 0 is considered 

{l/ri[Z] + a[B]) [j] = (l/? 7 [e*] + a[h^]) (6) 

called the CFIE, where g = y/ejl and generally a = 1/2 or 
ct = 1 in the literature (e.g. j^). 
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III. CFIE WITH HIERARCHICAL BASES 

The EFIE © produces ill-conditioned matrix [Z] when 
applied at low frequencies or on fine discretizations. This 
is due to the presence of the divergence operator and the 
opposite scaling of the two parts in A possible cure 

to this problem is based on Calderon’s identities ||7|. Another 
options is to use a hierarchical change of basis, with functions 
with increasing support. As described in Q, 0, this can be 
accomplished with a pre- and post-multiplication of the matrix 
[Z], obtained in standard RWG functions, by a matrix [Q] and 
its transpose, with a complexity scaling as N log N, where N 
is the number of unknowns. We refer to the hierarchical basis 
as Multi Resolution (MR) basis. This preconditioning effect 
is extremely efficient for very dense meshes, and it becomes 
less effective when the size of the support of the coarser 
functions is comparable with A/8, A being the wavelength 
of the surrounding medium. As described in fTQ| , the method 
is applied to large multiscale problems by: a) stopping the 
hierarchical construction at the desired level, and b) applying 
an algebraic preconditioner—as the LU or Incomplete LU 
(ILU) factorization—only on the part corresponding to the 
coarsest functions; that is precisely the regime where algebraic 
preconditioners are effective in smooth structures. The use of 
hierarchical loops as well as of hierarchical non-solenoidal 
functions is necessary to stop the hierarchical construction at 
a given level, and generate generalized RWG functions 
at the chosen coarse level; this, and the ability to efficiently 
find and use (possible) global loops constitute the dominant 
reasons for the use of hierarchical loops. 

From an algorithmic point of view, the MR scheme can be 
directly applied to the CFIE. It amounts to a pre- and post¬ 
multiplication of the Galerkin matrix in ([^ and a subsequent 
algebraic preconditioner, oblivious of the considered operator. 

However, the effectiveness of this approach is far from ob¬ 
vious. In fact, the identity operator appearing in Q questions 
the utility of using the MR basis with the CFIE. Indeed, while 
the MR change of basis effectively preconditions the EFIE, 
it could have a negative effect on the spectral properties of 
the operator of the MFIE, as explained in 0 0 From this 
perspective, the first question is if there is a trade-off between 
the regularizing effect on the EFIE part and the possible 
deterioration of the MFIE spectrum. 

As our results will show, in all investigated cases, including 
electrically large and complex structures, the use of the pro¬ 
posed hybrid hierarchical approach is always advantageous. 
We remark that the spectral properties of the CFIE are not 
simply the sum of the spectral properties of the constituent 
EFIE and MEIE; we also remark that the CEIE, in practice, 
only exists at the level of (finite) matrices and coefficients. 
A full theoretical explanation of this findings is beyond the 
scopes of this communication; the spectral properties of a 
similar form of the CEIE are studied, inter alia, in [?]; we 
are not aware of any other similar work on this topic. 

IV. Numerical examples 

In the first place, we verify that the MR basis does not 
corrupt the resonance free properties of the CEIE, considering 


frequency [MHz] 

EEIE 

CFIE 

RWG 

MR 

RWG 

MR 

300 

250 

275 

76 

57 

400 

278 

303 

72 

57 

474.56 

466 

470 

70 

60 

500 

250 

408 

70 

56 

600 

250 

540 

69 

56 


TABLE I 

Number of GMRES iterations to reach a residual lower than 

FOR THE RESONANT SPHERE TEST CASE, SECTIOnItvI 


an example with a sphere of radius 0.5 m illuminated by 
a linearly polarized plane wave at different frequencies. In 
particular, we hit the resonant frequency of 474.56 MHz. 
Table displays the number of GMRES iterations needed 
to reach a relative residual lower than 10“^. The CEIE 
formulation does not suffer around the resonant frequency and 
the MR change of basis preserves this property (as it was 
reasonable to be expected, as it does not alter the constituent 
equations) and actually improves the convergence rate. 

A. Role of hierarchical loops 

We focus on the preconditioning effect of the MR basis on 
complex structures with heterogeneous mesh and fine details; 
as reported in Section [T| an algebraic preconditioner on the 
coarsest level complements the MR preconditioning. Inter 
alia, this implies use of hierarchical loops. We start with the 
analysis of the dense-discretization breakdown, as a function 
of the edge refinement. We consider a cube of edge 1 m, 
at the frequency of 1 MHz, and meshes with decreasing edge 
length, denoted with h in Table We investigate the condition 
number and the number of GMRES iterations to reach a 
residual lower than 10“^, for an impinging linearly polarized 
plane wave. Table reports the results for the EEIE and the 
CEIE, both in the RWG, in the MR basis built with standard 
(point) loops, and in the MR basis built with hierarchical loops 
0. Consistent with results in 0, there is a slight improvement 
with hierarchical loops for the EEIE (with removal of a single 
very small eigenvalue). We note instead that the effect of 
hierarchical loops is significantly larger for the CEIE than for 
the EEIE. 


B. Role of the MR change of basis to handle details 

We apply the MR change of basis and the preconditioning 
strategy devised in fT0| to a closed object with a heterogeneous 
mesh with fine details. We consider a sphere of radius 0.5 m 
at the frequency of 300 MHz. Starting from a regular mesh 
with edge length of approximately A/10, we refine the mesh 
on part of the structure to have edges as short as A/200. We 
investigate the convergence as the level of detail increases. 
In Table |I^ we compare the number of GMRES iterations 
to reach a residual lower than 10“^, when the excitation is a 
linearly polarized plane wave. We use the standard RWG basis 
(RWG) and the MR basis with an algebraic preconditioner on 
the coarsest level (MR), both for the EEIE and the CEIE. The 
maximum (minimum) area of the mesh cells, (A™^^), 
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EFIE 



CFIE 


Ijh 

RWG 

MR hierarchical loops 

MR point loops 

RWG 

MR hierarchical loops 

MR point loops 

I 

3.40 10^ (14) 

5.5110° (9) 

1.12101 (9( 

1.25102 (jQ 

5.25 10° (8) 

1.39101 (10) 

2 

1.42 10^ (34) 

1.88101 (j5) 

5.86101 (i4( 

3.08 102 (33) 

1.20101 (j5( 

7.55101 (18) 

4 

5.23 10^ (51) 

3.92101 (21) 

3.37102 (2j) 

7.43 102 (54) 

2.24101 (20) 

4.31 102 (33) 

8 

2.07 10® (74) 

8.09 lOi (27) 

1.74 10® (30) 

1.66 10® (84) 

4.16101 (25) 

2.47 10® (70) 

16 

8.2910® (107) 

1.93 102 (37) 

6.09 10® (43) 

3.5610® (129) 

7.34101 (35) 

1.30 10^ (126) 

20 

1.39 10'^ (74) 

2.25 102 (41) 

1.94 10-1 (48) 

4.1010® (119) 

6.54101 (33) 

1.80 10^ (157) 


TABLE II 

Condition number (and GMRES iterations to reach a residual lower than 10“^) eor the cube example, Section IiV-AI 


^max y^min 

femax/(jmin 

EFIE 

CFIE 

RWG 

MR 

RWG 

MR 

9.4410° 

4.23 10® 

93 

21 

25 

21 

2.36102 

1.98 lOi 

186 

39 

45 

28 

1.0010° 

3.75 lOi 

212 

53 

59 

35 


TABLE III 

Number oe GMRES iterations to reach a residual lower than 

10“^ EOR THE EXAMPLE OE THE HETEROGENEOUS MESH, SECTIOn IiV-BI 


and the maximum (minimum) length of the mesh cells edges, 
^max Qf increasing heterogeneity of 

the mesh. Even in this simple case, we can appreciate the 
better convergence of the CFIE coupled with the MR and the 
algebraic preconditioner. 

V. Real-life cases 

Realistic cases are simulated with a MoM code implement¬ 
ing the MultiLevel East Multipole Algorithm (e.g. Q) to 
evaluate the Galerkin testing (ADE-EMS fTT|). We compare 
the results obtained with an iterative solver in the RWG basis 
with a simple diagonal preconditioner (RWG), in the RWG 
basis with a full near-field LU preconditioner (RWG LU), and 
in the MR basis with a LU preconditioner on the near-field 
interactions relevant to the coarsest functions |TQ| (MR), both 
enforcing the EFIE and the CFIE on the whole structure. 


A. ESA satellite test case at 2 GHz 

We consider a realistic test case by the European Space 
Agency, with a helix antenna on top (see Fig ure where 
mesh and size refer to the 12 GHz case of Section [W^. At the 
frequency of 2 GHz, the overall dimension is more than 130A. 
The mesh has edges from A/6 to A/130 long, corresponding to 
more than 10^ unknowns. When enforcing the CFIE, the EFIE 
is actually enforced on the open parts and on the wire antenna, 
which amount to nearly 0.2% of the entire functions. The 
iterative solver used is the BiCGStab, with a threshold on the 
residual norm of 10“^. Table 
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upper part, summarizes the 
results for different preconditioners and different formulations. 
The full history of the iterative solver is in Figure left panel. 
In this case, resorting to the CFIE formulation is mandatory 
to reach convergence, and the MR basis with an algebraic 
preconditioner on the coarsest functions saves more than half 
of the memory used in the RWG case, with a negligible 
increase in the number of iterations needed. 


B. ESA satellite test case at 12 GHz 

The same structure of Section [V^ is considered at 12 GHz, 
illuminated with the Spherical Wave Expansion of the feeder 
of one of the reflectors. The overall size is around 800A. The 
mesh has edges from A/4 to A/114 long, giving more than 
32 10^ unknowns. As the EFIE could not reach convergence 
in the 2 GHz case of Section V-A we test the CFIE only. In 
this very large case, the use MR preconditioning proves the 
only option to achieve practical convergence, as apparent from 
Table there, we see that we cannot reach convergence in 
the RWG basis and, on the other hand, the memory required 
to apply an algebraic preconditioner in the RWG basis is 
estimated to be more than 650 GB. Using the MR basis, 
we can apply the algebraic preconditioner only on the set of 
coarsest functions and keep the memory requirements feasible 
with the hardware available, which has 512 GB of RAM. The 
obtained surface current is depicted in Figure 


C. Low profile UAV 

We consider an unmanned aerial vehicle (UAV) with 
wingspan and tip-to-tail dimension around 30 A, illuminated 
with a linearly polarized plane wave. The structure is closed, 
with two deep cavities around the turbine. The considered 
mesh has edges from A/7 to A/65 long, resulting in around 
710+^ unknowns. The iterative solver used is the GMRES, 
with a threshold on the residual norm of 10“^. Table |lV| lower 
part summarizes the results for different preconditioners and 
different formulations. The full history of the iterative solver 
is in Figure right panel. The induced surface current is 
depicted in Figure In this case, convergence can be achieved 
with the EFIE formulation as well, with more iterations than 
with the CFIE. In the MR case, we reach convergence with 
some more iterations than in the RWG case with a full near¬ 
field LU preconditioner, but using less than half the memory. 
This applies to both the EFIE and the CFIE, meaning that the 
MR change of basis maintains the good convergence properties 
of the CFIE and allows an efficient algebraic preconditioner 
to be applied on a reduced part of the matrix. 


VI. Conclusion 

A hierarchical change of basis, developed to improve the 
convergence of the EFIE for dense discretizations, can be 
applied in conjunction with an algebraic preconditioner to 
increase the convergence of the CFIE for electrically large and 
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CFIE 


RWG 

RWG LU 

MR 

RWG 

RWG LU 

MR 

Satellite at 2 GHz Section 

V-A 








Ram [GB] 







13.1 

58.3 

24.0 

16.2 

62.3 

26.5 

BiCGStab iterations (threshold 

10-3) 

NC 

NC 

NC 

NC 

14 

27 

Satellite 12 GHz Section 

V-B 









Ram [GB] 







— 

— 

— 

285 

> 650 (extimated) 

477 

BiCGStab iterations (threshold 5 10 

— 

— 

— 

NC 

— 

99 

UAV Section 

V-C 












Ram [GB] 







19.4 

53.2 

25.6 

21.5 

48.1 

20.1 

GMRES iterations ( threshold 

10-3) 

NC 

538 

564 

NC 

71 

86 


TABLE IV 

Cases of Section|v] ‘NC’ means that no convergence is achieved in 1000 iterations, ‘—’ means the case has not been considered. 
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Fig. 1. Surf ace current, size, and mesh details for the satellite at 12 GHz of 
Section ES 




Fig. 2. Left panel: BiCGStab residual history for the case of S ection [V^ 
Right panel: GMRES residual history for the case of Section |V-C| Note 
that performances should be compared also in terms of RAM memory 
consumption, see Table [Tv] 
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multiscale structures. Leveraging on the hieararchical change 
of basis, the preconditioner can be applied on a reduced num¬ 
ber of unknowns maintaining the good convergence properties 
and the immunity to resonance problems of the CFIE. 

We have shown the advantage of the proposed approach on 
proof-of-principle examples as well as on complex, real-life 
simulations. 
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Fig. 3. Surface current for the UAV of Section [V^ with see-through slice 
and zooms on the blades of the turbine. 














































